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S£ ■ ABSTRACT 



We suggest that the growth of the central black hole in active galactic nuclei (AGN) 
due to the matter accreted over the AGN lifetime causes an evolution of the luminosity 
and spectrum. In a previous paper we have shown that the effective temperature of 
the UV continuum spectrum and the emission-lineequivalent width are anti-correlated 
with the black- hole mass. Here we estimate the change in the equivalent width of the 
£SJ ' emission lines due to the growth of the black hole mass caused by the accreted material 

o ' 
m 

t^. ■ models the evolving equivalent width is anti-correlated with continuum luminosity, 

implying that intrinsic evolution could contribute to the Baldwin effect. 

Q_i Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — quasars: 

Q ■ general — black hole physics — accretion disks 
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1. INTRODUCTION 

Quasars and Active Galactic Nuclei (AGN) are most probably powered by accretion of 



matter onto a supermassive black hole. The amount of accreted matter required to explain the 
the observed energy output of Seyfert galaxies and quasars is 0.001-10 solar masses per year, 
while the central mass required to maintain steady accretion for the observed luminosity (the 
Eddington limit) is of the order of 10 6 — 1O 9 M . Comparing these two numbers gives a timescale 
M/M ~ 10 7 - 10 9 years, which indicates a timescale for the evolution of the black hole mass and 
related properties, such as the emergent spectrum. 

As the black hole mass grows due to the accreted material, the observed luminosity and 
spectrum change with time. Since the continuum radiation ionizes the line-emitting material and 
drives the broad AGN emission-lines, the change in the properties of the radiation from the central 
source changes also the properties of the emission lines. In particular the spectral shape of the 
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ionizing radiation would determine the equivalent width (the ratio between the energy flux in the 
line to the continuum energy flux at the line frequency) of the emission-lines. 

The accretion disk postulated in AGN provides a convenient (although model-dependent) 
framework for which the emergent continuum spectrum may be calculated (e.g. Sakura and 
Sunyaev, 1973), depending on the black hole mass and the accretion rate. Alternatively, an 
empirical average AGN UV spectrum may be related to the characteristic black body temperature 
of the accretion disk . 

Baldwin (1977) has shown that the equivalent width of the CIV emission-line decreases with 
increasing continuum luminosity, a correlation known as the 'Baldwin effect'. A similar relation 
has been found also for other broad emission-lines. The origin of this effect is not understood (for 
a recent review see Shields and Osmer, 1998). 

In the first paper in this series (Wandel 1999a, hereafter paper I) we have argued that 
the Baldwin effect is related to the black hole mass. We have shown that the line equivalent 
width decreases with increasing mass, so that the expected correlation between black hole mass 
and luminosity in luminous AGN implies an inverse correlation between equivalent width and 
luminosity. 

In this work we suggest that at least part of the Baldwin effect can be attributed to 
evolution. We estimate the relation between the equivalent width and the continuum spectrum, 
and calculating the evolution of the continuum emission from the accretion disk for changing black 
hole mass and accretion ratewe show that for plausible evolutionary tracks the equivalent width is 
anti-correlated with the continuum luminosity. 

Section 2 gives an empirical relation between the model parameters (M and M) and the 
implied evolutionary time scales, as well as the observed spectrum and luminosity. In section 3 we 
suggest some plausible evolutionary scenarios for the accretion rate, and section in 4 we compose 
all these elements to produce the Baldwin effect . 

2. Accretion parameters and time scales 

2.1. The Evolutionary timescale for accreting Black Holes 

In order to maintain steady spherical accretion the luminosity must be less than the 
Eddington luminosity, 

L < L Edd = 4-7rGMm p c/o-T = 1.3 x 10 46 M 8 erg/sec 

or 

M 8 = 0.7^-^46 (1) 
where M 8 = M/1O 8 M , rj = L/L E dd is the Eddington ratio, and L 46 = L/10 46 erg/sec. 
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We may define the Eddington time t E as the time required for an accretion-powered object 
radiating at the Eddington luminosity to double its mass, due to accreted matter, 

t E = Mc 2 /L E = 4 x 10 8 y. 

The observed luminosity implies an accretion rate of M = 0.16e _1 L^qMqy -1 where e is the 
efficiency. Combining these two expressions gives the mass, accumulated if the accretion rate is 
maintained during a time t E '■ 

M 8 = 0.6e- 1 (L/L Edd )- 1 L 46 . (2) 
The timescale for evolution (for example, significant change in the mass) is given by 

t ev = M/M = €rf x T E 

for 7] = 0.1 this gives 3 x 10 8 yr for a Schwartzshild black hole and 1.6 x 10 9 yr for a Kerr black 
hole . 



2.2. The Mass-Spectrum Relation 

The continuum accretion disk spectrum depends on the black hole mass. The spectrum of the 
black body (and modified black body ) disk is given by integrating over the entire disk, 

[■Rout 

L v « / 2-KRB v [T(R)\dR (3) 
JRt 

where B U (T) is the Planck function and Rt is the transition radius from the optically thick to 
optically thin inner disk. 

In paper I we have demonstrated that the basic property of the accretion disk spectrum - 
thermal emission in the UV spectrum - leads to a peak in the SED, at an energy that is weakly 
decreasing with increasing black hole mass; From the approximate shape of the optically thick 
accretion disk spectrum 

vV 3 ( v 



L u « A — ) exp[ (4) 

(where A is a normalization constant and v co is the cutoff frequency) it is possible to derive an 
expression for the peak (or cutoff) energy of the UV spectrum: 

E co w E M 8 1/2 L% A (5) 

where E Q «17eV for an accretion disk around a Schwartzshild black hole and E a «20eV for a disk 
around a maximally rotating Kerr black hole , and E a is related to the maximal temperature in 
the optically thick part of the disk by E w 3kT max for a Schwartzshild black hole and 3.5kT for a 
Kerr one (Malkan 1990). 
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2.3. Relating the Black Hole Mass and Accretion Rate to the UV spectrum 

Plotting the cutoff frequency versus the luminosity it is possible to infer the mass and the 
accretion rate(or Eddington ratio; fig. 1 in Paper I). The L/LEdd ratio found in paper I is 0.01-1 
for quasars, 0.001-0.1 for Seyferts, while NLSls are all in the L/LEdd ~ 0.1 — 1 range. 

A similar analyses for the Schwartzshild accretion disk has been performed by Wandel and 
Petrosian (1988), who calculated for each pair of accretion disk parameters M and M the actual 
observables, e.g. the UV luminosity and spectral index rather than fitting the blue-bump. Wandel 
and Petrosian have calculated the accretion disk flux and spectral index at 1450A, for a grid of 
accretion disk parameters (M,M). Inverting the grid they have obtained contours of constant 
black hole mass and constant Eddington ratio (m) in the L — ajjy plane (fig 2 ) . Plotting in this 
plane samples of AGN it is possible to read off the contours the corresponding black hole mass 
and accretion rate. Comparing several groups of AGN a systematic trend appears: higher redshift 
and more luminous objects tend to have larger black hole masses and luminosities closer to the 
Eddington limit (see table 1). Similar results were obtained by Sun and Malkan (1989). 

Table 1: Grouping of AGN accretion disk parameters. 



AGN group 



LogF 145 o Log M/M & M/M Edd 



Seyfert galaxies 
Low Z quasars 



29-30.5 



28-29.5 



7.5-8.5 



8-9 



0.02-0.1 



0.01-0.5 



Medium Z quasars 
High Z quasars 



30-31.5 
31-32 



8- 9.5 

9- 9.5 



0.1-0.5 
0.03-2 
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Figure 1. Accretion-disk parameters and spectral properties in the a — L plane (Wandel and 
Petrosian 1988). Crosses: Seyfert galaxies, triangles, open circles and filled circles: low, medium 
and high redshift quasars, respectively. Continuous lines: constant mass; dashed lines: constant 
M/M (note that the labeling of the curves in the figure follows the notation rh = 16.7M / 'MecLcl)- 



2.4. The Equivalent Width - Mass relation 

In Paper I we have shown that the line luminosity for some lines - in particular CIV and 
OVI - is approximately proportional to the ionizing luminosity, as the line emissivity is weakly 
dependent on the ionization parameter. The luminosity in the line is given by 

Lline ~ LionfeffJO, (6) 

where f e ff is the effective covering factor, and j Q is the peak emissivity. The ionizing luminosity 
is determined by the maximal temperature of the optically thick part of the accretion disk , 
which has been shown to depend on the black hole mass approximately as M -1 / 4 for a variety of 
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accretion disk models. Combining these two results with eq. || above gives the mass dependence 
of the equivalent width - 

EW oc ea;j>[-(0.8 - 1.6)M g 1/4 ] (7) 

and similar dependencies for the luminosity (EW oc ex^CL^g 4 ). For the empirical power law 
AGN spectrum in the EUV (A < 1000A) we find EW oc with 7 ~ 0.13 - 0.25. 

3. Evolutionary Scenarios 

The evolution of the black hole mass over cosmological times depends on the accretion rate, 
which may of course change over time. From the theoretical efficiency and observed luminosity we 
can estimate the present accretion rate. In order to determine theblack hole evolution, we need to 
specify the variation of the accretion rate over time. We consider three representative scenarios: 

• constant accretion rate 

• Eddington limited accretion (constant Eddington ratio) 

• spherical accretion from a homogeneous medium 

The tracks in the E co vs. L plane of these three accretion scenarios are shown in fig. 2. See also 
the discussion in Wandel (1999b). 

3.1. Constant accretion rate 

This would correspond to external feeding of the black hole , which is regulated by an external 
potential, e.g. by stellar encounters, where the stellar motions are governed by the gravitational 
potential in the bulge of the host galaxy, far from the influence of the black hole . As the mass of 
the black hole grows, M does not change, and hence the accretion parameter, m ~ M/M ~ M . 
Eventually accretion rate will become very sub-Eddington , and the black hole will be starved, or 
at least on a diet. Since L ~ M, in this scenario the (bolometric) luminosity also is also constant. 
Note that the spectral energy distribution does change (though slowly) as M increases. 

3.2. Eddington limited accretion 

This would be the case in a black hole that is over-fed. The accretion rate cannot become 
very super-Eddington , and the source will regulate the accretion rate to L i=a LEdd orma e -1 . 
Since ra ~ M/M, m ~ constant implies M ~ M and hence L ~ M. 
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3.3. Spherical accretion from a homogeneous medium 

In this case we assume the matter supply comes from a homogeneous distribution (gas or 
stars), due to the gravitational potential of the black hole . On large scales, the accretion will 
be spherical, and the accretion radius at which the black hole gravitational potential becomes 
significant is given by 

R acc « GMvl « 3M 8 vl 00 pc, (8) 
where v* = 300t>3ookm/s is the stellar velocity dispersion. 
The accretion rate is given by the Bondi formula 

M « 47ri4^* 2 P* = (O.3M /yr)M| U3 - o 3 o ( 10M ^ pc _ 3 ) > ( 9 ) 

where is the stellar number density. 

In this case L ~ M ~ M 2 and m ~ M, that is, the Eddington ratio increases with 
time. Eventually the Eddington ratio will approach unity and the accretion will become 
Eddington limited. 



3.4. Accretion from an inhomogeneous medium 

A similar expression can be derived for a non-homogeneous medium, with a density profile 
~ R~ p . Assuming the velocity dispersion is independent of R, eqs. || and |9| give for the 
accretion rate in that case 

M ~ M 2 - p . 

If also the velocity dispersion has a radial dependence, v* ~ R~ q then we get a more complicated 
dependence, 

M ~ M^~ q ~ p . (10) 

In this case, any functional dependence is possible. If for example (7 = 1/4 and p = 1/2, 
M ~ M° , and if g = 1/2 (point mass) and p = 1/2, M =const. 
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Figure 2. Black Hole evolution for the three accretion ratescenarios described in the text in the 
cutoff energy vs. the monochromatic luminosity at 5100A. Diagonal dashed lines indicate constant 
black hole masses, and diagonal dotted lines indicate the Eddington ratio. Vertical arrows describe 
the direction of evolution for a constant accretion rate, Horizontal arrows show the evolution for 
a spherically accreting black hole from a homogeneous medium, and diagonal dashed arrow are 
for accretion at a fixed Eddington ratio. The superimposed data is PG quasars (diamonds; from 
Boroson and Green, 1994), Seyfert 1 galaxies (triangles and squares, the latter indicate AGN with 
BLR reverberation data from Wandel, Peterson and Malkan 1999, with luminosities modified for 
H = 50) and Narrow Line Seyfert Is (Boiler, Brandt and Fink, 1996) 



4. Putting All Together: an Evolutionary Baldwin effect ? 



In the previous section we have found the relation between the accretion parameters M 
and M. It is easy to find the functional dependence of the continuum luminosity on the black 
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hole mass for several accretion scenarios. This cab be shown as evolutionary curves in the M — rh 
and UV-spectrum vs. luminosity plane (figs. 1 and 2). For example, an Eddington bound 
accretion disk will evolve along the constant rh curves (dashed curves in fig. 1). A spherically 
accreting black hole will evolve along curves of rh ~ M, which correspond to nearly horizontal 
(from left to right) lines in fig. 1, and a constant accretion rate will yield rh ~ M _1 which are 
almost vertical (upward) lines in fig. 1. Similar evolutionary tracks are shown in fig. 2 (note 
however that the vertical axes of figs. 1 and 2 have an opposite sense). 

Combining this with the dependence of the effective radiation-temperature on the accretion 
parameters, and using the result that the emission-line luminosity is positively correlated with 
the radiation temperature, it is possible to find the relation between the equivalent width and the 
continuum luminosity for each pair of accretion scenario and temperature model. 

4.1. The L-T relation 

The dependence of the effective temperature on the black hole mass is summarized in table 
2. The header shows the three accretion scenarios and the corresponding dependence of the 
continuum luminosity on the mass. The next row gives the dependence of the luminosity on the 
mass for each scenario. The left column shows the temperature models. For each pair the table 
gives the functional dependence of T e ff on M. 

Since the black body -radiation boundary and the photospheric models, have a very similar 
dependence on the accretion parameters, we combine them in one row in tables 2 and 3. Table 
2 has nine (T e ff,M) pairs, for accretion disk models, and an additional row for the model 
independent black body - variability estimate of T e ff (marked "T(BB-var)"). 

For all pairs except two (constant M with T(BB — rad) and homogeneous medium accretion 
with Tbb)- The effective disk temperature is decreasing with increasing black hole mass. We note 
that this is true also for the general temperature-luminosity relation (eqs. 21,22 in Paper I), so 
the accretion disk model is not essential. In that case the anti-correlation between temperature 
and luminosity is built into the the expression for the temperature. 

For all accretion scenarios the optical and UV continuum luminosity is increasing with M. 
This is certainly the case for Eddington bounded accretion and for accretion from a homogeneous 
medium; For a constant accretion rate the bolometric luminosity is constant, but as the black 
hole mass increases the peak of the accretion disk spectrum moves to lower frequencies. Since 
the optical or near UV continuum are on the Rayleigh-Jeans part of the black body spectrum, 
the luminosity in these bands increases as the mass increases even if the total luminosity remains 
constant. However, since for a radially extended accretion disk the low energy slope in the near 
UV is much flatter than for a Planck spectrum, this is a very small effect. As the M=const is a 
limiting case, and in reality the accretion rate is expected to grow, even if slower than linearly, 
as argued above, we replace the M=const. column in table 3 by a more general case, 
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f3 < 1, which allows for an evolution in luminosity. 

Applying the dependence of L on M it is possible to deduce the dependence of T on L. 
Relating the EW to the spectral UV peak energy (Paper I) we can quantify the expected Baldwin 
effect for each of the models. 

Decreasing accretion rate. In the three cases considered above, the accretion rate increases 
with time or remains constant. What if the accretion rate decreases with time? In the two latter 
accretion disk -temperature models in Paper I (5.2.c and d) we find E max ~ M" 1 / 3 , so T e ff 
increases, and so does the line luminosity. If we take L oc M the continuum luminosity decreases, 
and there is still a negative correlation. This is the case also in the general L-T relation. In the 
model, the temperature is independent of the accretion rate, and since the mass increases, the For 
the black body and radiation-pressure dominated models the situation is less clear, because T e ff 
decreases with decreasing M (as M 1 / 4 /M 1 / 2 for the black body model and being independent of 
the accretion ratein the radiation-pressure dominated model, it still decreases as M -1 / 4 ). Since 
L ~ M decreases, the question whether T e ff is correlated or anti correlated with L would depend 
on the details of the behavior of the accretion rate. 



How does the equivalent width vary? The equivalent width is defined as Lu ne /F\, where F\ 
is the flux per unit wavelength at the wavelength of the line. In Paper I we have shown that Lu n , 
depends on the ionizing continuum, which correlates strongly with T e ff. The luminosity in the 
line is also roughly linearly dependent on the continuum luminosity, L ~ \F\, so that 



f(T) being a monotonic function with an exponential increase near 3kT ~ lRyd. This gives for 
the equivalent width 



Applying the observed UV spectral slope, L v oc v~ a (1.5 < a < 2), we can derive from 
the Temperature - Mass dependence shown in table 2 a relation between equivalent width and 
luminosity. From Paper I we have 



so we can derive the Baldwin effect expected for each model quantitatively. Table 3 assumes 
a = 2. For a L v ~ v~ a (a < 2) spectrum the exponent of L should be multiplied by a — 1. 

We see from table 3 that for most combinations of effective temperature model and accretion 
scenario there is an evolutionary Baldwin effect , and the predicted slope (dlnEW / dlnL) is 
7 = 0.25(a — 1), consistent with the observed Baldwin effect of 0.15-0.2. 



4.2. The equivalent width 



Lline — f(T e ff)L, 



EW(L)~L line /L = f[T(M)}. 
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4.3. The Exceptions 

We conclude that for many accretion-scenarios and temperature-models the line luminosity 
decreases with increasing mass, except the two for which the effective temperature is constant, 
implying an evolutionary Baldwin effect : as the black hole evolves over cosmological times due to 
the mass accreted, the continuum luminosity (optical or near UV) increases, while the equivalent 
width decreases. Also the two pairs with a constant effective temperature will eventually show a 
Baldwin effect . To see this, note that for the (Tbb — M(spherical)) pair, m is increasing as M, 
so it will eventually approach the Eddington accretion rate, and move from the Tbb model to one 
of the other T e ff models, which (for M (spherical)) do have a Baldwin effect . In the other pair, 
(constant M with T(BB — rad)), the situation is opposite, but with the same result: since there 
m ~ M _1 , eventually the accretion rate will become enough sub-Eddington so that the Tbb 
model will apply, which for the constant accretion rate scenario does have a Baldwin effect . 

5. discussion 

In paper I we have made the case for a Baldwin effect due to the dependence of the equivalent 
width on the black hole mass. We have shown that the line equivalent width decreases with 
increasing mass, so that the expected correlation between black hole mass and luminosity in 
luminous AGN implies a Baldwin effect . 

In this work we show that an intrinsic evolution of the black hole mass and accretion ratecan 
create an evolutionary Baldwin effect . In addition, the relation between the emission-line 
luminosity and the ionizing continuum shape depends on the unknown density profile of the line 
emitting gas, and on the unknown differential covering factor, which could show intrinsic evolution 
(for an individual object along it's lifetime) as well as cosmological evolution. 

How can we infer whether evolution is the main cause of the Baldwin effect? This work 
shows that evolution is a possible cause. But there are several other possible causes as well. For 
example, the ranges in actual black hole masses and accretion rates in individual AGN might 
deviate greatly (and systematically) from the average relations assumed. 

6. Summary 

We suggest that the growth of the central black hole mass due to the accreted matter 
causes the luminosity and the spectrum to evolve over over the AGN life time. Using general 
model-independent relations, or the thin accretion disk spectrum, we estimate the evolution of the 
equivalent width and the continuum spectrum, and show that for plausible evolutionary tracks as 
well as for the model-independent black body temperature estimate and for most variants of the 
thin accretion disk model the equivalent width decreases with increasing continuum luminosity 
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implying an evolutionary origin to the Baldwin effect. 

Stimulating discussions with Jack Baldwin during the meeting "Quasars as Standard Candles 
in Cosmology" at La Serena, Chile, triggered this work. I also acknowledge valuable discussions 
with Gary Ferland, Matt Malkan and Vahe Petrosian, the contribution of an anonymous referee, 
and the hospitality of the Astronomy Department at UCLA and Stanford University. 
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Table 2: The relation between effective accretion disk temperature and black hole mass for various 
accretion scenarios 



Accretion scenario: Constant M Constant L/LEdd Spherical 



L(M) oc M const ~ M ~ M 2 
T e ff model: 

T BB M- 1 / 2 M-V4 C onst 

T rad M' 1 / 4 M' 1 / 4 M" 1 / 4 

T(r* = 1) or T(BB/rad) const M' 1 ' 4 M~ 4 / 2 

T(BB-var)cx L-V4 con st M^ 1 / 4 M" 1 / 2 
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Table 3: The predicted Baldwin effect (dependence of the EW on L) for different effective 
accretion disk temperature models and accretion scenarios. An EUV spectrum of L v ~ v~ 2 has 
been assumed. For a L u ~ v~ a (a ^ 2) spectrum the exponent of L should be multiplied by a — 1. 
m_i = m/0.1. 

Accretion scenario: L ~ M@ L/L 

Edd Spherical 

m : m ~ const ~ M 

______ ______ 

Tge 4mi / 1 2 L7 6 1/4 L (2/3-3)/4 M (i-/?)/2 L -i/4 congt 

T rad 4mi / 1 4 L7 6 V4 L (/3-2)/4 M (l-/3)/4 £-1/4 £-1/8 

T(r* = l) a 3L 4 " 1/4 L-V4 L -i/4 L -i/2 

T(BB-var) b lOL^ 4 L -1 / 4 L -1 / 4 L^ 1 / 4 



"and the T(BB/rad) boundary model 



